Passing quantum correlations to qubits using any two-mode state 
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We draw an explicit connection between the statistical properties of an entangled two-mode continuous 
variable (CV) resource and the amount of entanglement that can be dynamically transferred to a pair of non- 
interacting two-level systems. More specifically, we rigorously reformulate entanglement transfer process by 
making use of covariance matrix formalism. When the resource state is Gaussian, our method makes the ap- 
proach to the transfer of quantum correlations much more flexible than in previously considered schemes and 
allows the straightforward inclusion of the effects of noise affecting the CV system. Moreover, the proposed 
method reveals that the use of de-Gaussified two-mode states is almost never advantageous for transferring en- 
tanglement with respect to the full Gaussian picture, despite the entanglement in the non-Gaussian resource 
can be much larger than in its Gaussian counterpart. We can thus conclude that the entanglement-transfer map 
overthrows the "ordering" relations valid at the level of CV resource states. 

PACS numbers: 03.67.-a, 03.67.Mn, 42.50.Dv, 03.67.Lx 
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The establishment of a long-haul entangled channel is one 
of the most important and challenging aims of that part of 
quantum technology concerned with communication and dis- 
tributed computation [1]. The amount of experimental and 
theoretical efforts that have been put into the design and real- 
ization of a device allowing such a result with a sufficient de- 
gree of reliability is tremendous |2[] . Among the rich plethora 
of schemes put forward so far, an intellectually stimulating 
and pragmatically inspiring proposal relies on the possibil- 
ity to quantum mechanically correlate two remote systems by 
transferring somepre-available entanglement initially carried 
by a resource Jj, 0, |5[] . This stems as a rather natural way of 
achieving a distributed channel for quantum communication, 
especially when the entangled resource is embodied by the 
continuous variable (CV) state of a photonic system, which 
has been proven to be a reliable courier of quantum correla- 
tions. In fact, such an idea has been implemented, recently, in 
a setup involving cold atomic ensembles and light fields 
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Despite the scheme is attracting an increasing attention also 
in virtue of its versatile nature, which makes it suitable for im- 
plementations in various settings, including cavity and circuit 
quantum-electrodynamics, the formalism itself is largely in 
need of development. In its basic form, the protocol requires 
the arrangement of a bilocal interaction between each remote 
qubit and one of the modes of the CV photonic resource. The 
aim is to arrange a situation suitable to the pouring of quan- 
tum correlations from the CV resource to the non-interacting 
two-level systems in such a way that, upon tracing out the two 
modes, we are left with an entangled state of the qubits. In this 
procedure, it is often the case that the explicit form of the CV 
resource, decomposed in some basis, is needed. This makes 
the computational efforts required for the evaluation of the de- 
gree of transferred entanglement quite demanding. And yet, it 
is rather sensible to believe that asking for the full CV state de- 
composition is probably not at all necessary, especially if one 
considers the class of Gaussian entangled two-mode states, 
which are fully determined by assigning only a few of their 
statistical properties. This is exactly the main point of the 



present investigation: we develop a self-consistent formal ap- 
proach that is able to reconstruct the entries of the two-qubit 
density matrix resulting from the entanglement transfer pro- 
cess simply by means of statistical properties of the pre-built 
resource. This represents the first step of a theoretical inves- 
tigation directed towards the generalization and extension of 
the panorama of situations where entanglement transfer can 
be conveniently and almost effortlessly studied. This includes 
also the experimentally important case of non-Gaussian CV 
resources having the form of s-photon subtracted two mode 
states 1211. The proposed formalism allows us to quantitatively 
study the performance of this class of states in relation to the 
tasks discussed here. 

The remainder of this paper is organized as follows. In 
Sec. |T] we lay down the notation and recall some basic defi- 
nitions and tools useful for the treatment of CV systems. In 
Sec.|n]we attack the case of entanglement transfer from a re- 
source prepared in a general Gaussian state. We give explicit 
expressions for the elements of the two-qubit reduced density 
matrix resulting from the transfer process as functions of the 
statistical moments of the resource. We highlight their theo- 
retical and computational convenience by addressing a gen- 
eral two-mode squeezed thermal state, which allows us to un- 
veil the spoiling effect, on the degree of transferred entan- 
glement, of the resource's thermal character. A similar study 
is then performed with respect to a dissipation-affected two- 
mode squeezed state: in both cases we demonstrate that the 
degree of poured entanglement is set once the resource is as- 
signed. Sec. [Ill] tackles the case of experimentally relevant 
non-Gaussian states obtained by subtracting s photons sym- 
metrically from a two-mode squeezed state 121 la] . A series 
of results are achieved in this case: first, we demonstrate that 
the Gaussian core of the photon-subtracted resource (i.e. its 
statistical moments up to the second) completely specifies the 
amount of transferred entanglement. However, the use of such 
a non-Gaussian state is not always advantageous with respect 
to its Gaussian counterpart (obtained by not subtracting pho- 
tons): despite photon-subtraction works as an entanglement 
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distillation protocol for the CV state, a two-mode Gaussian 
state is almost always a more efficient resource for entangle- 
ment transfer. We also draw some conjectures on the relation 
between the discrepancy in the performances of Gaussian- 
associated and original non-Gaussian states and the degree of 
non-Gaussianity of a resource H]. In Sec.[[V]we summarize 
our findings and highlight the procedure for adapting our tech- 
nique to the recently proposed protocol for entanglement re- 
ciprocation [10]. Finally, in Appendices I and II we provide 
a few technical steps needed in order to fully understand our 
formal approach. 

L PRELIMINARIES 

We consider a bipartite CV system comprising modes j = 
1, 2, each described by its respective bosonic annihilation and 
creation operator cij and fit, such that [fij, a' k ] = Sjk with Sjk 
the Kronecker symbol being 1 for j = k and otherwise. For 
the aims of this work, it is worth introducing the quadrature 
operators Xj — ctj + fit and yj = i(fit — fij), which are 

collected in the vector of field quadratures Q = (x\ yi x^yi)- 
A generic two-mode state with density matrix pyi can be 
completely described in phase space by its Weyl characteristic 
function iflTll 

X(£,fj) = Tr[p 12 ^! (0^2 fa)], (1) 

where £ = £ r + i& and rj = r\ r + ir/i are the phase-space vari- 
ables, and we have introduced the phase-space displacement 
operator Dj(a) — exp[aa| — a*a,j] {a E C). It is straightfor- 
ward to show that pi2 is explicitly related to the characteristic 
function as 11211 

P12 = 4, / d 2 Zd 2 r)x(t,r))Di(-0®D 2 (-V)- (2) 

The exact characterization of a generic CV state requires in 
principle the knowledge of all the infinitely-many moments of 
the quadrature operators. The first moments can be adjusted 
by local displacements and are thus totally irrelevant for the 
purposes of characterizing and processing entanglement (they 
are assumed to be set to zero throughout this paper without 
loss of generality). On the other hand, second quadrature mo- 
ments are often crucial: they can be conveniently collected in 
the covariance matrix V which, for a two-mode system, has 
entries = Tr[pi 2 {Qi,Qj}/2]. The covariance matrix of 
any physical two-mode state can be cast, via local unitary op- 
erations (which do not affect entanglement by definition), into 
the standard form 1 1311 

V = ( ni m ) (3) 
\m n 2 J 

where each of the two block matrices nj = Diag[nj, rij] ac- 
counts for the statistical properties of the corresponding mode 
j (jij > 1), and m = Diag[m + , m_] specifies the cross-mode 
correlations (m+ > > 0). From now on, unless other- 

wise specified, we assume that the covariance matrix of any 
state at hand has already been put in standard form. 



A very special role in the CV arena is played by Gaussian 
states. Gaussian states are defined as having a Gaussian char- 
acteristic function ITill 

X(t,v)=e-^ V T , (4) 

with Q = (£ r r\ r rji) the vector of canonical phase-space 
variables. Therefore, an arbitrary two-mode Gaussian state 
P12 is completely specified (up to the first moments) by as- 
signing the covariance matrix V. The density matrix is di- 
rectly expressed as a function of the covariance matrix via 
Eqs. (f2]l and @J. This greatly simplifies the mathematical 
treatment of Gaussian states, which despite living in a infinite- 
dimensional Hilbert space are effectively described by a finite 
number of degrees of freedom. 

Physically, bipartite entanglement in a generic CV state 
is "signalled" by the violation of the Heisenberg uncertainty 
principle by the partially transposed density matrix Il4 llal . 
For arbitrary non-Gaussian states, such an inseparability con- 
dition can be cast in terms of a hierarchy of inequalities in- 
volving moments of arbitrary order jlol . In the special case 
of Gaussian states, clearly, entanglement is completely quali- 
fied and quantified by algebraic combinations of the elements 
of the covariance matrix V only. Specifically, a two-mode 
Gaussian state is entangled if and only if the following in- 
equality is fulfilled, 

v- = -^=\/a- VA 2 -4detV < 1 (5) 
v2 

with A = dctrii + dctri2 — 2dctm iTfl [nil . In general, 
Eq. <(5j stands as a sufficient condition for entanglement in an 
arbitrary non-Gaussian state with covariance matrix V. 

H. BIPARTITE GAUSSIAN-STATE CASE 

In this Section we will derive a complete formal description 
of the entanglement transfer process from two-mode Gaussian 




FIG. 1: (Color online) Scheme of principle for an entanglement 
transfer process. An entangled two-mode CV state is generated 
with an off-line process. Each mode addresses a two-level system 
qj (j = 1,2). Through bilocal interactions (represented by the cir- 
cular arrows in the figure), part of the entanglement initially shared 
by modes 1 and 2 is transferred to the joint state of systems q\ and 
q2- The CV modes can then be discarded or can enter a detection 
stage for outcome post-selection, like in the entanglement reciproca- 
tion scheme IIOTI . 
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states to two-qubit systems in terms of the covariance matrix 
of the Gaussian resource. Eq. (f2j) is the starting point of the 
present analysis: we now use twice the decomposition of the 
identity in terms of Fock states X)n° m =o l n > TO )i2 ( n > m \ = 1 
and utilize it in order to get the Fock-state decomposition 
P12 = E^m=oE^g=o7^|n,TO)i2(p,5| with the coeffi- 
cients 



^pi = — 

him _2 



(6) 



Following Refs. JH d 



\&'/*c^\\Z\ 2 ) 



(7) 



where / np (£) = x (ra|£>i(-£)|p)i. 
flSll . one has that for n > p 

fn P (0 -- 

with (x) an associated Laguerre polynomial of degree s 
and argument x (r 6 Z). Analogous explicit forms hold for 
frnqiv)- Although the four-fold summation involved in the 
full decomposition of pi2 seems daunting, we now show that, 
once the interaction model between the remote qubits and the 
two-mode CV system is assigned, the entanglement transfer 
scheme allows for some simplifications that considerably re- 
duce the computational difficulties. 

For the sake of definiteness, we assume that the two-level 
systems are both prepared in their fundamental energy level, 
associated with the logical state \g) qj (j = 1,2). Any other 
choice is, obviously, equally valid. Then, in order to remain 
along the lines of the research conducted so far on entangle- 
ment transfer, we consider bilocal Hamiltonian models of the 
form 



n 



(.7 = 1,2) 



(8) 



with the x, y-Pauli matrix a^ v and the coupling rate SI, which 
we take to be the same for both qubit-mode subsystems (again, 
an assumption that can be easily relaxed). Eq. (0 embod- 
ies a resonant dipole-like coupling Hamiltonian under the 
standard rotating wave approximation, an interaction mech- 
anism of broad relevance that has attracted a considerable 
body of theoretical and experimental work in the last fifty 
years [19j]. We stress that the validity of the arguments that 
will be presented throughout this paper are not bound to 
this specific choice. Any coupling Hamiltonian could well 
be taken, mutatis mutandis, without affecting the main im- 
plications of the present study. The full form of the time- 
propagator e~ lHj ' q j t , decomposed in the single-qubit basis 
{|e) g ., \g) qj }, has been given in 12011 . Here, it is sufficient 
to state that when exactly n photons populate mode j, one 
has \g,n) gjtj -> C n (T)\g,n) gjtj - ig w (r)|e, n - l) qj j with 

r = Sit, C n {r) = cos(t^) and S n (r) = y/1 - C*(t). 
Equipped with these tools, one can now track the effective 
evolution of the two qubits obtained upon partial trace over 
the CV degrees of freedom. More formally 

p qiq , (r) = Tr 12 [e- 1 ^ T P i 2 ®\gg) qiq2 (gg\ e ^ T ] . 

(9) 



A simple calculation allows us to determine the 7^ coef- 
ficients that are associated with each entry of the two-qubit 
density matrix p qiq2 {r). One finds the relevant correspon- 
dences presented in Table [I] whose form can be verified upon 
explicit calculation (the Hermitian conjugates entries are eas- 
ily found). The entries of Table [fl are the only coefficients 
that are not identically null. This leaves us with the following 
two-qubit density matrix 

/Anmir) G nm (T)\ 

B nm (r) V nm (r) 

„tlo ° Vnm{T) Cnm{T) ° 

(r) £ nm (r)J 

(10) 

with £ nm {r) = 1 - Anm(r) - S„ m (r) - C nm (r). The form 
of P\2{t) is significant. The coherence term J2 n m ^nm( T ) 
arises in virtue of symmetric processes where both the qubits 
simultaneously absorb or emit an excitation, being thus re- 



sponsible of the transformation \gg) q 



ee 



On 



the other hand, T> nm (r) accounts for anti-symmetric events 
where while q\ emits or absorbs an excitation, under- 
goes the opposite physical process, thus giving rise to the 
lff e )<3i?2 l e 3)<?i<?2 transition. Such parity-related effects are 
reflected into the apices and pedices of the associated 7^ co- 
efficients, whose form is still to be determined. This is the 
task of the following analysis. 

It is important to stress that 7^'s are numbers depending 
merely on the statistical properties of the two-mode CV sys- 
tems via the Weyl characteristic function x(^v)- As such, 
they can be written as explicit functions of the elements of the 
covariance matrix V associated with a given Gaussian CV 
state. Here, we highlight the relevant technical steps required 
for such calculations. In order to discuss a significant exam- 
ple, we consider the case of 7™™. First, one can conveniently 
change variables as £ = re 1 ^ , -q = se %e (with r, s S [0, 00 [ and 
4>, 6 G [0, 2ir]). The four-fold integral entering the definition 
of 7^* is thus changed into a double integral 



nm 
/ nrn 



-1 poo />oo 

^ / rdr / S d S e-^' 2+s M V)4°V) 



(i 







x d<p d9x{re i4, ,se w ). 
Jo Jo 

(11) 

The double integral involving the angular variables can 
be worked out by writing explicitly the form of the 
characteristic function and using a power-series expan- 



TABLE I: Coefficients entering p qiq2 (t), Not reported are the coef- 
ficients that turn out to be identically null upon explicit evaluation. 
For simplicity of notation, we have dropped any r-dependence. 



Density-matrix entry 




Coefficient 


Iss)<7i92 (gg\ 




s~i'2 s~i'2 nm 
m — ^n^mlnm 


\9e) qiQ2 (ge\ 


B 


f~i2 o2 nm+1 
— u n J m+l7nm+l 


l e ff>9i<?2< e 5l 




q2 ^~»2 n + lm 
— °n+l^m7n+lm 


\gg)ni2{ ee \ 


Qnm ~ — Cn 


C C Q _.n + lm+l 
,&n-\-l^rn,&rn-\-l Jnm 


\ge) qiq2 (eg\ 


T^n m — Cn Sn + 1 Cm Sm-\-\ 'Y n m _|_ \ 
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sion. We end up with d(j)J^ d9 x( re ^ , sel9 ) = 
4^e- " ir2 ^ 2s2 EZo( rs ) 2k+1 Gk(m±) with 

. (2fc-l)!! , fc /l , m 2 -m^ 

Gfe(TO±) - ki»! m 2Fi U- 15 ) - 

(12) 

where 2 Fi is the Gauss hypergeometric function. The re- 
maining radial part of the integral is easily evaluated by 
exploiting the power-series decomposition of the associ- 
ated Laguerre polynomials 112 111 and the Gaussian integral 

/ °° dre'^r 2k+2i+1 = (k + i)\/2 Mi € Z JH. After some 
simplifications, one gets the final expression 

°° 4m 2fc / 2 \ 

^ [(ni + l)(n 2 + l)f +1 V "1 + 1/ 

/ 2 \ / 1 to 2 - mi \ 
x 2 Fi -m, 1 + fc; 1; — 2 F X -, -fc; 1; . 

(13) 

Any other 7^ involved in p gig2 (r) can be calculated follow- 
ing similar steps, which we omit here. In Appendix A we give 
the explicit form of all the relevant ones. It is worth stress- 
ing that, regardless of the apparently intimidating form of the 
coefficients in Eqs. $13[ and Eqs. (IA-U - (IA-3b . they are rather 
useful and pragmatic. Their plug-&-play nature allows an un- 
precedented flexibility in the study of entanglement transfer. 
It is enough to assign the covariance matrix of the two-mode 
Gaussian resource, without the necessity of knowing the full 
decomposition of the state onto a given basis, in order to un- 
ambiguously determine the amount of entanglement that can 
be transferred to two remote qubits prepared in their ground 
state. In the next Sections we demonstrate the flexible na- 
ture of our results by addressing significant quantitative ex- 
amples. It is usually the case that the infinite sums involved 
in the determination of p qiq2 (r) can be truncated using appro- 
priate cutoffs, their value depending on the features of V. For 
instance, for a two-mode squeezed state of squeezing factor 
(elwe have m,2 = cosh.(2£),m± = ±sinh(2£). Obvi- 
ously, Tri 2 (pi 2 ) = 1, which translates into J2 n ,m inm = 
as it is easy to check. In order to explicitly evaluate the nor- 
malization of the state, we have truncated the sum over k in 
Eq. ( fT3l and analogous at k c — 100 while those over n and m 
at n c — m c — 25. These cutoffs provide excellent results for 
values of ( up to ~ 1.5, as seen in Fig.|2](a). Larger squeezing 
parameters require larger n c and m c values, as it is intuitive to 
understand. It is well-known, in fact, that the probability that 
the Fock state with n photons is populated in a squeezed state 
of £ > 1 is non-negligible for quite a large range of n values. 

The validity of the results presented so far can be conve- 
niently tested by using this very same instance of two-mode 
CV resource for entanglement transfer and comparing the out- 
comes with available literature on the subject yfl. Through- 
out this work, we use the entanglement measure based on the 
negativity of partial transposition criterion lfl5ll . which is de- 
fined as twice the modulus of the negative eigenvalue of the 
partially-transposed two-qubit density matrix with respect to 
one of the two qubits M23I1 . Using the coefficients calculated 
here and setting ( = 0.86, which is known to optimize the en- 



tanglement transfer with this type of resource iH, we find the 
behavior shown in Fig.[2](b), which is in perfect quantitative 
agreement with the results of Refs. tfH] and thus states the reli- 
ability of our findings. It should be stressed that, as discussed 
in Ref. the aperiodic behavior shown in Fig. (f2]l (b) re- 
sults from the interference of pure Rabi oscillations, each oc- 
curring at frequency y/nfl, associated with transitions within 
each two-level system induced by the chosen multi -photon re- 
source. The incommensurability of the periods of such oscil- 
lations is responsible for the destruction of any periodicity. 



A. Thermal nature and dissipation effects in the CV resource 

We now start discussing the remarkable flexibility of the 
established connection between transferred entanglement and 
statistical properties of the CV resource. Our first step is the 
use of a more general Gaussian resource than a standard two- 
mode squeezed state. The most general state of a single-mode 
CV system is given by a displaced squeezed thermal state. 
However, when the multi-mode scenario is considered, many 
possibilities of combining and ordering single-mode as well as 
two-mode unitary operations are available. In order to include 
a standard two-mode squeezed state as a limiting case, here 
we concentrate on the following situation: we consider single- 
mode squeezing (along arbitrary directions in phase space) 
of two CV modes in thermal equilibrium at their respective 
temperature. Such squeezed thermal modes are then superim- 
posed at a beam splitter of transmittivity T. More formally, 
we consider the state 

Pl2 = B(Ri ® i? 2 )(5i ® &)pS(& ® ^(Rx ® &)t# 

(14) 

with B = cxp[i cos _1 (v / T)(a;i?; 2 — yix 2 )] the beam splitter 
operator, Sj and Rj the phase-space single-mode squeezing 
and rotation operators 111 ill (whose form is here omitted). We 
have introduced the tensor product of two thermal states — 
5^ m=Q Pnm \n,m)i2{n,m\, which we take to have mean 
thermal occupation number n% and n 2 respectively, so that 



(a) (b) 




FIG. 2: (Color online) (a): Normalization of a two-mode squeezed 
state evaluated as X]" c Em" 7nm> shown against the squeezing pa- 
rameter f . The calculations are performed at fixed cutoffs k c — 100 
and Tie = m c — 25. (b): Entanglement passed to two remote 
qubits upon bilocal interaction with a pure two-mode squeezed state 
of squeezing factor ( — 0.86, plotted against the dimensionless time 
r = Qt £ [0, 2tt]. The calculation has been performed retaining the 
cutoffs given in panel (a). The result perfectly matches the studies in 
Refs. (1. 
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Pnm = + 1)" +1 (^2 + l) m+1 [11]. Depending 

on the beam-splitter transmittivity and the relative direction of 
squeezing, the mixed state pi 2 can show some entanglement. 
While its Fock-state decomposition is somehow uncomfort- 
able and definitely lengthy to compute, it is immediate to find 
the form of the covariance matrix V. Each of the unitaries in 
Eq. (fl4l > has, indeed, a symplectic counterpart that can be used 
to transform the covariance matrix of the initial thermal state, 
which reads V p th = (2ni + 1)1.2 © (2n 2 + 1)1 2 with l k the 
k x k identity matrix. One has Sj = Diag[e~% e 

cos (fij 1 2 + i sin (fij<Ty and a = 



Rj — 

-vT^Tl 2 

Here, sj and ipj are the single-mode squeezing parameter and 
phase-space rotation angle, respectively. The covariance ma- 
trix thus changes as 

V P12 = B T {R 1 ®R 2 ) T (S 1 ®S 2 ) T V pth (S 1 ®S 2 ){R 1 ®R 2 )B. 

(15) 

Eq. j!5\ , which is readily calculated upon simple matrix prod- 
ucts, is then put in standard form by following the general 
recipe given in Ref. 01311 . With this, the evaluation of the trans- 
ferred entanglement for a given choice of the parameters en- 
tering Eq. ( fT3T > is immediate. In order to study the effects that 
mixedness in the CV resource has in the transfer performance, 
a step which has only been superficially addressed so far, here 
we discuss the case of si = —s 2 = 0.86, tp 12 — 0, T = 1/2 
and ni,2 = n. Interesting effects arising from the use of 
asymmetric two-mode Gaussian resources (i.e. H\ ^ n 2 ) 
will be addressed elsewhere UM . With these choices, for 
n = we have the covariance matrix of a two-mode squeezed 
state and we thus retrieve the results of Fig. [2] (b). At non- 
zero thermal occupation number, on the other hand, we have 
«i,2 = (l + 2n) cosh( 1.72) and m± = ±(l + 2n) sinh(1.72), 
which lead to Fig. [3] The amount of transferred entanglement 
decreases with n, as it should be expected given that the entan- 
glement within the CV resource itself is spoiled by an increas- 
ing mixedness n ^ 0. However, the relevant feature is that 




6^0.0 



FIG. 3: (Color online) Transferred entanglement against dimension- 
less interaction time r and the mean thermal occupation number 
ni,2 = n of the two CV modes used as resources. For the whole 
range of n shown in this plot, the initial CV resouce is nonseparable. 



we find large temporal regions of separability of the qubits' 
state, despite the CV resource remains entangled within the 
full range of values of n shown in the plot. In fact, we get 
2/_ = 1 (cfr. Eq. (O) at n = (e 2 ^ — l)/2, which is equal to 
2.292 for ( = 0.86. The thermal nature of the CV state thus 
makes the entanglement-transfer process unreliable. 

On the other hand, one can also consider the effect that 
dissipation affecting the CV resource has on the passage of 
quantum correlations. Again, having related the transferred 
entanglement to the elements of the covariance matrix of the 
two-mode state proves to be a major advantage. In fact, by 
assuming weak coupling between each mode j = 1,2 and its 
respective bosonic bath at thermal equilibrium, which allows 
our analysis to be kept within the Born-Markov approxima- 
tion, the dissipative dynamics of the CV system is described 
by the master equation B25I1 

r 2 _ 

dtPu = - \N(2ajp 12 a,j - p 12 }) 

j=i (16) 

+ {N +\){2a ]Pl2 a] - {ojoj.pu}) 

with r _1 the lifetime of a photon in the dissipative environ- 
ment characterized by its mean thermal number N. By mak- 
ing proper use of standard operator correspondences 112511 . this 
master equation can be changed into a dynamical equation for 
the characteristic function v)- After a tedious but other- 
wise straightforward calculation, one gets 

F 2 

dtxit V) = ~ g E^ 1 + 2N ^ I' + + Mi%M£> V) 

(17) 

with jiti = £ and \i 2 = r), Eq. (flTT i is readily solved and, by 
means of definition @, it is finally found that the initial co- 
variance matrix V(0) evolves towards the covariance matrix 
of the bath (2N + 1)1 4 as 



V(t) = (2JV + 1)(1 



e- rt )ll4 



V(0)e 



-n 



(18) 



which is in agreement with the analysis reported in Refs. [27 
We now plug the covariance matrix of a two-mode squeezed 
state as V(0) in Eq. ( [T81 and study the modifications induced 
in the entanglement transfer function, for set values of N and 
£, by an increasing "dissipation time" Ft. This dimension- 
less parameter is the product of the dissipation rate V, which 
characterizes the environmental action, and the time t dissi- 
pation has acted on the CV resource before the entanglement 
transfer process begins 02611 . The results are shown in Fig. |4] 
where £ = 0.86 is taken with N = 0.1. Quite expectedly, the 
transferred entanglement decreases with Ft and disappears for 
rt ~ 0.5. This matches what is found for the two-mode CV 
resource, whose v_ becomes larger than 1 for 



Tt > log 



'47V(7V + 1) + sinh(2r) - cosh(2r) + 1 



47V 2 + 47V 



(19) 

which equals 0.52 for the choice of parameters made above. 
One can thus claim that the entanglement transfer process is 
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FIG. 4: (Color online) Qubit entanglement transferred from a two- 
mode squeezed state subject to dissipation, plotted against the di- 
mensionless interaction time r and dissipation time Ft for N = 0.1 
and C = 0.86. 



effective as far as the CV resource is entangled, although re- 
gions of separability appear, dynamically, within the temporal 
window studied in Fig. [4] This behavior has been checked 
to hold true also when the initial two-mode covariance ma- 
trix describes a general physically-allowed state other than a 
two-mode squeezed vacuum. 



HI. BIPARTITE NON-GAUSSIAN RESOURCES 

The handy nature of the connection established between 
statistical properties of the entangled resource and transferred 
two-qubit entanglement turns out to be even more striking 
when non-Gaussian CV states are at hand. In this Section, we 
show that the initial Gaussian component of a de-Gaussified 
two-mode state is key in the determination of the entangle- 
ment shared by the remote two-level systems. We consider 
de-Gaussification as obtained by photon-subtraction from a 
general Gaussian CV state lf7 H28ll . 

An intuition of our claim comes from considering the def- 
inition of an s-photon subtracted two-mode CV state, which 
reads 

p\2 = ^ J ^Ax(^)(a?£i(-0^*)(^-^)4*) 

(20) 

with M the normalization factor and s the number of photons 
subtracted (symmetrically) from the two modes of the CV re- 
source. Although this equation is only formal, it encompasses 
the crucial features entailed by the physical de-Gaussification 
process of subtracting a photon. In Appendix B we also high- 
light the results achieved upon use of an effective procedure 
for subtracting photons, making use of highly-biased beam- 
splitters and photo-resolving detectors. 

It is straightforward to check that the determination of the 
coefficients associated with the various elements of the two- 
qubit density matrix resulting from the interaction with the 



non-Gaussian CV resource tracks the steps depicted in Sec. [TT] 
and Appendix A. In fact, all one has to do is to take 

pq (n + s)\(m + s)\(p + s)l(q + s)\ +sq+s 

lnm V n\m\p\q\ In+sm+s K > 

in the entries of Tab. U and matrix p qiq2 . The pres- 
ence of the square-root factor reflects the effects of de- 
Gaussification process. This result demonstrates our claim 
that the entanglement-transfer capabilities of an s-photon sub- 
tracted state, under the model for local Hamiltonian addressed 
within the context of this work, are determined by the knowl- 
edge of the Gaussian "core" of the state and, of course, the 
number of photons being subtracted. 

Our task now becomes manyfold. First, we aim at show- 
ing that a non-Gaussian resource obtained from a two-mode 
squeezed state, which is a realistic and interesting case to 
study [7], beats the corresponding Gaussian resource only 
for proper choices of the squeezing parameter £. However, 
the maximum value of entanglement transferred by a de- 
Gaussified state for r 6 [0, 2tt], which is a reasonable time 
to wait, never exceeds the one achieved via a Gaussian re- 
source. This strongly contrasts with the entanglement of the 
state itself: the negativity of an s-photon subtracted two-mode 
squeezed state is larger than the one for s = 0, regardless of 
( > and s > 0. In the remainder of this Section, we indi- 
cate with £ Gauss (SnonGaus) the transferred negativity when 
a Gaussian (de-Gaussified) CV resource is used. In Fig.|5]we 
show the difference between the entanglement £cuas trans- 
ferred upon usage of a Gaussian two-mode squeezed state and 
ZnonGaus when a 1-photon subtracted two-mode squeezed 
state is employed. Such a difference is evaluated at set val- 
ues of the dimensionless interaction time r and the degree of 
squeezing of the Gaussian-core resource. We notice a rich 
structure of maxima and minima. Notably, with the exception 




FIG. 5: (Color online) Density plot of the difference Eghus — 
SnonGaus between the entanglement transferred using a two-mode 
squeezed state and that achieved via a 1-photon subtracted two-mode 
squeezed state. The transferred entanglement is studied against the 
dimensionless interaction time r and the squeezing parameter £. The 
side color-bar indicates the color scale. 
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FIG. 6: (Color online) Density plot of the difference Egcius — 
SnonGaus between the entanglement transferred using a two-mode 
squeezed state and that achieved via a 7-photon subtracted two-mode 
squeezed state. The transferred entanglement is studied against the 
dimensionless interaction time r and the squeezing parameter £. The 
side color-bar indicates the color scale. 



of the three dark areas in the region of small £, witnessing 
that the use of the non-Gaussian resource favors the transfer 
of entanglement to the remote two-level systems, such a dif- 
ference remains largely non-negative throughout the entire re- 
gion determined by r G [0, 2tt] and £ S [0,2]. The results 
are thus clear: the use of a Gaussian CV resource is (almost 
always) much more effective than a non-Gaussian state (be- 
longing to the class of de-Gaussified states studied here) when 
entanglement transfer processes are addressed. This qualita- 
tive behavior holds for any value of s tested by our calcu- 
lations. For instance, for s = 7 we get the plot in Fig. [6] 
which shows that the three negative areas where the use of a 
de-Gaussified state is more advantageous get squashed in the 
low-squeezing region, while the moduli of Scans - ZnonGaus 
increase, thus showing an even more striking convenience 
in using Gaussian states. The additional square-root factor 
in Eq. ( l2T1 i. which deeply affects the interferences among 
components associated with fixed numbers of excitations, is 
the reason for such a considerable change in behavior of the 
discrete-variable entanglement function. The squashing to- 
wards regions of lower squeezing originates from the modifi- 
cations induced by photon-subtraction on the photon-number 
probability distributions of the states being used: as s grows, 
such probability distributions are peaked at smaller £ and their 
width shrinks with respect to what occurs for the s = case, 
thus explaining the behavior observed in Figs. [5] and [6] In or- 
der to show that Gaussian states allow for the achievement of 
the largest transferred entanglement, regardless of the form of 
the entangled resource, we have used the general formalism 
of Sec. Ill Al to generate a 1000-element sample of bona fide 
random covariance matrices. They have been used in order to 
evaluate the maximum qubit entanglement (for r € [0, 2tt]) 
that can be achieved by using such fully Gaussian resources 
and the s-photon subtracted states having the latter as Gaus- 
sian core (with s = 1, ..,4). The result is that, quite clearly, 



Gaussian states are able to achieve the largest discrete-variable 
entanglement, as far as the model for bilocal interaction dis- 
cussed here is considered. Moreover, the successive subtrac- 
tion of photons reduces the maximum of entanglement being 
transferrable. These features are clearly shown, for a subset 
of only 22 elements, in Fig. [7] where the transferred entangle- 
ment upon use of one of such elements is studied against the 
number of photon subtractions performed in the CV resource: 
the decrease with s, regardless of the Gaussian core part of the 
resource, is quite evident. 

We now address a second interesting point related to the use 
of non-Gaussian CV states. General considerations in the the- 
ory of CV entanglement reveal that the Gaussian state having, 
as covariance matrix entries, the moments (up to the second) 
of a general two-mode CV state, provides a lower bound to 
the entanglement content of the latter. In this sense, Gaussian 
states are said to be "extremal" lT29ll . Recently, it has been 
proven that the entanglement in an s-photon subtracted two- 
mode squeezed state is upper- and lower-bounded by a func- 
tional of the second moments of such non-Gaussian state l30ll . 
Here, we show that this is actually not the case for the amount 
of transferred entanglement: the one achieved via such "fic- 
titious" Gaussian resource (with equal first and second mo- 
ments) is often larger than the entanglement passed to two re- 
mote qubits by the (in general non-Gaussian) CV state. This 
can be very easily checked by proceeding as follows. One can 
take the covariance matrix elements of an s-photon subtracted 
two-mode squeezed state as given in Ref. 13011 , use the expres- 
sions valid for the Gaussian-resource case (cfr. Eq. (fT3l and 
the formulas in Appendix A) and compare the results with 
what is gathered by considering the formal apparatus, pre- 
sented earlier in this Section, for s photons being subtracted 
from it. The explicit calculations confirm the general trend 
anticipated here that Gaussian CV resources appear to be op- 
timal for entanglement transfer purposes. Moreover, it allows 
us to conclude that extremal states in the CV scenario are not 
mapped, in general, into extremal discrete-variable states by 




4 



FIG. 7: (Color online) We show the maximum amount of entangle- 
ment transferred within r g [0, 2tt] from a sample of 22 randomly 
generated s-photon subtracted Gaussian state for s = 0, ..,4. As 
more photons are subtracted from the Gaussian core part of the state, 
despite an increasing resource entanglement, the maximum quantum 
correlations that can be passed to two remote qubits decreases. 
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the entanglement transfer process at hand. 

We are currently looking for a physically clear relation- 
ship involving such a discrepancy between a non-Gaussian 
state and its Gaussian-equivalent one and the degree ofnon- 
Gaussianity of the starting C V resource as quantified by quan- 
tum relative entropy J9[]. For the case of an s-photon sub- 
tracted two-mode squeezed state, such a figure of merit sim- 
ply coincides with the von Neumann entropy of the Gaussian- 
equivalent state associated with it. We already have a few 
numerical evidences showing that the squeezing-dependent 
functional form of the difference between the maxima of 
transferred entanglement in the non-Gaussian and Gaussian- 
equivalent cases mimics the shape of the degree of non- 
Gaussianity in the low-squeezing region. 



IV. CONCLUSIONS AND PERSPECTIVES 



modes, labelled for simplicity 1 and 2, which are initially in a 
separable state. Upon proper projection of the qubits, an ex- 
pression analogous to ( 1221 would be obtained, with TL\$ to be 
interpreted, this time, as qubit projectors. Clearly, the effec- 
tiveness of the reciprocation of entanglement would depend 
on the sort of projection being chosen. We conclude that the 
results discussed in this paper can be fully exploited in order to 
assess the case of entanglement reciprocation as well, which 
depend crucially on the very same coefficients 7^ calculated 
in this work. 

We expect that the handiness of the results achieved through 
our methodology will trigger further development and deep- 
ening of entanglement-transfer processes that, in light of re- 
cent experimental progresses along these very same lines fl, 
hold the promises to embody a pragmatically viable route to 
long-haul distribution of quantum correlations. 



We have provided a general formalism for the treatment of 
entanglement-transfer processes from an arbitrary two-mode 
resource to remote two-level systems. The remarkable flex- 
ibility of the method proposed here allows for the general- 
ization and extension of this protocol for the distribution of 
long-haul quantum communication channels to situations that 
have been only partially addressed so far, such as the use of 
thermal/dissipated Gaussian CV resources and experimentally 
interesting de-Gaussified states 

The proposed methodology is quite versatile and can be 
straightforwardly adapted in a way so as to consider multi- 
mode CV resources, as in Ref. Jj], and include the effects 
of post-selection and detection, such as in the entanglement 
reciprocation scheme of Refs. fiol l3lll . In this context, the 
entanglement transfer protocol y|] as described and studied 
here is complemented by "hard" projective measurements per- 
formed on the entangled resource. One can consider two com- 
plementary cases. In the first, a two-mode CV state is used to 
entangle two separate qubits by means of local interactions 
and post-selective measurements of a specific and appropriate 
nature. While Ref. [1 1 Ofl considered the case of a projection 
of the CV resource onto coherent states of a proper ampli- 
tude m, 

other choices, such as parity or homodyne mea- 
surements, are possible (see Ref. 03211 for a recent example). 
In general, upon projection of modes 1 and 2 by means of the 
operator III II2 such that ITj = Ilj (J = 1, 2), it is imme- 
diate to see that the reduced density matrix of the two qubits 
is cast into the form 



Pq±q 2 



E pq n np 
I nm 



e"f( n i 



(22) 



1j t ) [p^ 9 (n2,T)] a time-dependent operator 



with 0™f(II 

spanning the Hilbert space of qubit qi ((72) whose form, for 
an assigned interaction Hamiltonian, depends on the choice 
of III (II2) and on the initial preparation of the qubits. The 
coefficients 7^ are calculated as described in the previous 
Sections. One can also consider the reverse situation where 
the state of qubits qi and q2, entangled as a result of the pro- 
cess described above, is used as a resource to entangle two 
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APPENDIX A 

In this Appendix we give the explicit form of the 7^ co- 
efficients involved in the two-qubit density matrix resulting 
from the interaction with a CV system prepared in a Gaussian 
state. In Eq. (fT3l we have already provided the one associ- 
ated with the I gg) qi q2 (gg | projector. We start looking at the 
coefficient for \ge) qiq2 (ge\, which depends on Jnm+i- By 
definition 



nm+l _ 

mm-fl 2 



d 2 ^d 2 rj x(i,V)fnn(Ofm+lrn+l(v) 



4m 



21; 



^ [(m + i)(n 2 + i)p+ 
/ 

2F1 -to- 1,1 + k; 1; 



r2 *i -n, 1 



fc;l; 



m + 1 



n 2 



1*1 



1 m_ 
— , k\ 1; 



(A-l) 

In a completely analogous way, we get that 7™+i„ is given by 
Eq. (lA-lb with the replacement n — > n + 1 and m + 1 — > m. 
As for the \gg) qiq2 (ee| element we have 



n+lm+l _ 
inm 



A k+1 [{k + l)!]V(m + l)(ra+l) 



K k (m±) 



x 2*1 -n,2 + fc;2 



(2fc+l)![(m + l)(n 2 + l)] fe + 2 

2 \ „ / _ . _ 2 



1 



2*1 -n,2 + fc;2; 



n 2 + 1 
(A-2) 

with Kk{m±) a combination of the hypergeometric func- 
tions 2*1 (±5> —k] 2; 1 — ~r-^j with coefficients entirely de- 
termined by m±. A somehow analogous form holds for the 
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coefficient associated with \ge) qiq2 (eg\, which reads 



„ +lm _ 2 2fc + 3 [(fc + l)!]V(m+l)(n + l) 



' nm+1 



(2fc + l)![(a +!)(&+!)] 



fc+2 



L k (m±) 



x 2 Fi ( _ n , 2 + fc; 2; ) 2 -Fi (-m,2 + fe;2; 



a + 1 



6+1 
(A-3) 

As before, Lk(m±) is fully determined by the correlation 
terms of the CV systems's covariance matrix and contains hy- 
pergeometric functions. 

All the remaining density matrix elements are associated 
with coefficients that turn out to be identically null. In or- 
der to give an intuition for this, let us consider 7," T ™ +1 , which 
accounts for \gg) qiq2 (ge\. The angular part of the four-fold 
integral to be performed reads 



account of the effects of this operator onto a Gaussian state 
and, then, onto the entanglement transfer process, here we as- 
sess the fully physical case. For the sake of simplicity, we 
address the single-photon case only. This is done by consid- 
ering 



P12 



A'' 



d 2 ^ 2 ^x(C^) 3 (i|B 13 i)i(-O|0) 3 



4 f 3 |i) 



x 4 (l\B 24 D 2 (- v )\0) 4 (0\Bl 4 \l) 4 , 



(B-l) 



where we can write Di{— £)|0)3(0| = 

£jr*=o/«p(OIM>ia<P,o| and £ 2 (-»?)|o>4<o| 

El, F o fmqiv) |m, 0) 2 4<<Z, 0|. One can then use the 
Fock-state decomposition of the state resulting from the 
action of a beam splitter ITTIl . For instance, we have that 



d8 X (re l4 \se l9 )e- l ° = 



(A-4) 



for any covariance matrix V in standard form as in Eq. (0. 
Similar arguments hold for any other identically-zero coeffi- 
cient in p qiq2 (t). 



APPENDIX B 

In this Appendix we show the formal derivation of the co- 
efficients valid for the case of a photon subtracted state, when 
the implementation of the non-Hermitian operator &j occurs 
via highly-biased beam splitters and the detection of single 
photons. We assume that the modes 1 and 2 of a bipartite 
Gaussian state are superimposed to two additional modes, 3 
and 4 respectively, at beam splitters of transmittivity T. For 
T — * 1, the beam splitters tap at most a single photon from 
modes 1 and 2, this occurring with a small (but non-zero) 
chance. The tapped photons populating modes 3 and 4 are 
then revealed at two single-photon resolving detectors. Upon 
measurement of one photon per mode, a single-photon sub- 
tracted state is achieved. By cascading s of these stages, an 
s-photon subtracted state is conditionally produced and the 
operator afa^ is implemented. While Eq. d20l > gives a formal 



B 13 \n,0) 13 = J2(-^ k TH^-T 2 



k=0 




By putting everything together, we eventually get 

(1-T 2 ) 2 



Pl2 



n,p—0 m,q—0 



X V(n + l)(m + +l)(q + l)f n+lp+1 (0fm+i q +i(v) 
_ (1 - T 2 ^ 2 



-^E 



oo oo 

m+1? 



n,p— m,q— 



(B-3) 

The factor (1 — T 2 ) 2 , which does not depend on the sum- 
mations' indices, is washed out by the normalization fac- 
tor N = tt 2 /[(1 - T 2 ) 2 E n , m (n + l)(m + lfrStf^i 1 ]. 
From this point on, the time dependence of the two-qubit den- 
sity matrix elements resulting from the entanglement trans- 
fer process can be evaluated as highlighted in the body of the 
manuscript. For T — > 1 any of them becomes identical to the 
corresponding expression valid for the formal case, as studied 
in Sec. [Till Quantitatively, for T > 99.99% the formal and 
physical approach give results that are indistinguishable from 
each other. 
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